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Abstract. We consider a random variable Y and approximations Y n , n € N, denned on the 
same probability space with values in the same measurable space as Y . We are interested in 
situations where the approximations Y n allow to define a Dirichlet form in the space L 2 (Py) 
where Py is the law of Y . Our approach consists in studying both biases and variances. 
The article attempts to propose a general theoretical framework. It is illustrated by several 
examples. 
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I. Bias operators. 



We consider a random variable Y defined on a probability space (Q, A, P) with 
values in a measurable space (E,^) and approximations Y n , n G N, also defined on 
(Q,A,f) with values in (E,^). In the whole study we suppose that there exist an 
algebra V of bounded functions from E into R or C dense in L 2 (E, J 7 , Py) containing 
the constants and a sequence (a„) n eN of positive numbers, with which the following 
hypotheses are considered : 



Vy? G V, there exists A[<p] G L 2 (E, T , Py) si. Vx G £> 
lim n ^ ooQ;n E[(^(y n ) - <p(Y)) X (Y)] = E Y \A[<p]x]. 



the expectation Ey being relative to the law Py. 

( VcpeV, there exists A[y>] G L 2 (£, .F,Py) si. Vx G X> 
(H2j i lim^oo a n E[(v9(F) - <p(Y n )) X (Y n )] = E Y [A[ip] X }. 



(H3) 



Vy? G X>, there exists j%] G L 2 (£, J 7 , Py) si. Vx GjD 

Um n ^ 0O a n E[( ¥ »(y n ) -^(y))( x (y„) -xQO)] = -2E y L%]x]. 



We first note that as soon as two of hypotheses (HI) (H2) (H3) are fulfilled (with the 
same algebra T> and the same sequence a n ), the third one follows thanks to the relation 

~ A + A 



When defined, the operator A which considers the asymptotic error from the point of 
view of the limit model, will be called the theoretical bias operator. 

The operator A which considers the asymptotic error from the point of view of the 
approximating model will be called the practical bias operator. 

Because of the property 

< A[<p], X > L 2(Py)=< <f, A[X\ >L2(Py) 

the operator A will be called the symmetric bias operator. 

Remark 1. Under (HI) the limit of a n E[((p(Y) — (p(Y n ))x(Y n )] exists and is equal to 

E Y [A[x]p]-E Y [A[p X ]]- 

The operator A with dense domain possesses an adjoint A . If V C T>(A ) then (H2) 
is satisfied and 

(l) A[p] = T[<p] - <pT[i\. Vpev 

Reciprocally, if (HI) and (H2) are supposed and if 1 G T>(A ), the map x l—> Ey [A[</?]x] = 
EyL4[x]<^] — Ey[A[v?x]] is continuous and so is x ^ E[]4[x]</?] which shows V C VA 
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and relation (1) holds. We see that the hypothesis 1 G T>(A*) is rather strong, it will 
be not fulfilled in general. 

Example 1.1. Let us take for (EjJ 7 ) a metrisable compact set with its Borel a-field 
and let (X t ) be a Feller process with values in E and transition semi-group (P t ) (A Feller 
process on an Led. space reduces to this situation by the Alexandrov compactification 
see [12] chap. XIII §20-21). 

Let (Ap,Dp) be the generator of the strongly continuous contraction semi-group 
(P t ) on C(E). Suppose (P t ) be in duality with a strongly continuous semi-group (Q t ) 
with respect to a probability measure v and let (Aq,Dq) be the generator of (Q t ) on 
C(E). Then, if there is an algebra of bounded functions containing constants V C 
D p DDq dense in L 2 (u), the approximation X t of X satisfies hypotheses (HI) to (H3) 
and we have on £>: _ 

A[(p]= A P [(p] 

A[<p] = Aq[<p] - <pAq[1] 
A*[<p] = A[<p]+<pAq[1]. 

Indeed, if tp, x G V 

jE u [(<p(X t ) - (p(X )) X (X )] = j<P tV - V , X >u -> <A P [(p],x >u 

^E„[(<p(X )-<p(X t ))x(X t )} = 1[< Q t <p-<p+<p(l-Q t l), X > v ] ^< A Q [ V ]- V A Q [llx>u ■ 

Hence (HI) to (H3) are fulfilled and, by theorem 1 below, the limit \~E v [((p(X t )— (p(X ) 2 ] 
extends to a symmetric Dirichlet form on L 2 (u). o 

The basis of our study is the following theorem : 

Theorem 1. Under hypothesis (H3) 

a) the limit 

QJ 

S[p, X ] = lim ^E[(v9(r n ) - tp(Y))(x(Y n ) - X (Y)] 

defines a closable positive bilinear form whose smallest closed extension is denoted 
(£,B). 

b) (S,B>) is a Dirichlet form 

c) (£, D) admits a square field operator T satisfying V</?, X £ ^ 

r[p] = A[<<?\ - 2<pA[<p] 
E Y [r[tp} X } = lima n E[(tp(Y n ) - V (Y)) 2 ( X (Y n ) + X (Y))/2] 

n 

d) (£, D) is local if and only if\/ip G V 

lima n E[( V (Y n )- V (Y)) 4 ] = 0. 

n 
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Demonstration, a) That (£, V) be closable comes from the Friedrich construction of 
the minimal selfadjoint extension of a symmetric operator. Let us recall the argument. 
By £[np, x] = — < A[(p],x >l 2 (p y ) Vy?, x e ^ the form (£, V) satisfies 

u n G V, \\u n \\ ^ =>- £[u n ,v] ^ Vt; G D 

and this property suffices to imply closability (cf. [13] ex. 1.1.2, [9] Chap. I ex. 1.3.4 
or [5] lemma 111.24). 

b) In order to prove that the form (£,B) is Dirichlet we will use the following 
elementary property : 

If K is a compact subset ofR, Ve > there exists a polynomial p(x) such that 

(i) < p(y) — p(x) < y — x Wx < y G K 

(ii) \p(x) -x\ <e Vr G [0, 1] n K 
(Hi) p(x) > — e Vx 6 K. 

Let (-Ra)a>o be the strongly continuous contraction resolvent associated with (£,D), 
we have to prove that the operators \R\ are sub-Markov (cf. [13], [9], [23]). For that, 
since here \R X 1 = 1, VA > 0, because £[l,tt] = V« G B, it is enough to show that 
< u < 1 R x u > and this for A > 1 since R a = ± J2T=i(( 3R *+p) k Va, /3 > by 
the resolvent equation. 

Let u be a measurable function from E into 1st. < u < 1, denoting as usual £\ 
the form £ + X\\.\\ 2 L2 , the functional 

11 

<b{v) = £[ v ] + X\\v - -|| 2 o6D 

A 

satisfies 

®(R x u) +£ x [R x u-v] = $(v) veB. 

It follows that R x u is the unique element of D minimizing <3> on D, and the sequence 
v n converges in © (equipped with the norm \fE{) to R x u if and only if lim n $(i> n ) = 
®(Rxu). 

Since V is dense in D, let v n G V be such a sequence converging to R x u. Let p n be 
a polynomial fulfilling property (i) to (Hi) for e = and X containing the range of the 
bounded function v n . Let us put w n = p n o v n . Since V is an algebra, w n G T> and by 
(i) we have £[iu„] < £[v n \. Now, by (%) and (ii) \p n o v n - y\ < \v n - y\ + ± Vy G [0, 1] 
so |w n -f| < |f„-f| + ^ what gives |K t -f|| 2 2 < \\v n - ^|| 2 2 + ^\\v n - f\\ L 2 + ± and 

j « j . 2 .. u.. 1 

$K < $ K + - K - t + 

n A rr 

The sequence t>„ being bounded in L 2 , it follows that $(u> n ) - ► ^K-Ra - ") and u>„ — > i?A« 
in D. Taking eventually a subsequence converging a.s. and using w n > — - by (in), we 
obtain .Raw > what proves the property. 

c) Let us denote (A, VA) the selfadjoint operator associated with {£, D) (Friedrich's 
extension of (A,T>)), since the algebra V C PA is dense in D, the theorem 4.2.2 



4 



of [9] applies and the form (£,D) admits a square field operator satisfying T[(p] = 
2[y? 2 ] - 2(pA[(p] V(/9 E T> and T[ip] = A{ip 2 ] - 2^A[ip] if tp EVA and ip 2 G PA The 
formula of the statement comes from 

E Y [T[<p]x] = E Y [A[<p 2 ] X ] - 2E Y [A[<p]tp X ] 

coming back to the definition of A. 

d) Let us remark that for f E V we have 

4£[f, f] - 3£[f 2 , f 2 ] = Yima n E[(f(Y n ) - f(Y))% 

n 

So, if the assumption of the statement holds, using the fact that one dimensional 
contractions are continuous on D (cf. []) 

4£[u 3 ,u] -3£[u 2 ,u 2 } = VKGDnl 00 

this is enough to imply that 8 is local (cf [9] Chap. I §5). Reciprocally, if £ is local, 
since £ admits a square field operator, the functional calculus applies (cf. [9] Chap. I 
§6) and 

A£[u 3 ,u] - 3£[u 2 ,u 2 } = 2E[3u 2 T[u]] - -E[Au 2 T[u]] = 
VueDn L°° and the condition of the statement is fulfilled. o 

Comment. Considering for Y a Brownian motion B indexed by [0, 1] as random 
variable with values in C([0, 1]) and taking for Y £ the approximation Y £ = B + ^feW 
where W is an independent standard Bronian motion, we may apply the theorem with 
V the linear combinations of functions ip(B) = e 1 ^ f dB with regular / say C\. 
We have with x(B) = e l ^ 9dB 

E[(e^o /dye - e^fo f dY )(e i -fo 9dY E _ e iQ g dY^ 

= E[e i ^ s+9)dY ]E[(e i ^i ' fdW - l)^ 1 ^^ gdW - 1)] 

so that 

\vra -EMY £ ) - V (Y)(x(Y £ ) - x(Y)} = (- f fg &) e -*/o</+»> a * 
£ Jo 

what may be written —2 < A[ip], x > with 

m=j""[-\J fdB-\J m 

as seen by an elementary calculation. Hypothesis (H3) is satisfied. The theorem yields 
the well known Ornstein-Uhlenbeck structure on the Wiener space (see e.g. [9] or [27]). 

We can say that from a pedagogical point of view, in order to introduce the error 
calculus on the Wiener space (basic Malliavin calculus) and the same would be true for 
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the Poisson space or the Monte Carlo space (cf. [5]), theorem 1 is a quite convenient 
tool. It is simpler than the theorem on infinite products of Dirichlet structures used in 
[9] or [5]. It allows also to construct Dirichlet forms in a variety of situations as will 
show the examples below. 

This being said, when studying mathematically approximations, the most usefull 
part of the theorem is often the easiest one to prove, i.e. part a), because the limit form 
is often recognized as a classical form whose properties (Dirichlet character, square field 
operator, locality) are known. 

Remark 2. Suppose that instead of (H3) we assume that the limit lim ra a n E[(ip(Y n ) — 
(p(Y) 2 ] exists Vtp G T>. Then, denoting 2S[ip] this limit, if the form (S,T>) is closable, 

i) the hypothesis (H3) is equivalent to V C VA where [A, VA) is the generator of 
the form (£,D) smallest closed extension of (£,£>) (Indeed, u G VA 44> / i— > £[/, u\ is 
continuous on V, hence if P C T>A hypothesis (H3) is satisfied with A = AonV and 
if (H3) is satified V C VA) 

ii) parts b) and c) of the demonstration of the theorem apply and show that (£, D) 
is Dirichlet with square field operator. 

Remark 3. Under (H3) the condition d) of the theorem Vy? G V lim n a n E[((p(Y n ) — 
tp(Y)) A } = is equivalent to either of the conditions : 

(j) 3A > 2 lim n a n E[|^(y n ) - <p(Y) \ x ] = \/tp G V. 

(jj) VA > 2 lim n a n EMY n )- V (Y)\ x ] = V</? G V. 
Indeed, it suffices to proves that (j) implies (jj). Let A be such that (j) is true, then for 
H > A 

a n E[\ V (Y n ) - tp{Y)\»\ < a n E[\ V (Y n ) - V 9(r)| A ](2|| V 9|| 00 )^ A - 
and for A > /j, = 2 + e 

a n E[\(p(Y n ) -^(y)H < (a n E\(p(Y n ) - ^(Y)\ 2 ) 1/2 (a n E\^(Y n ) - ^(Y)\ 2+2e ) 1/2 

iterating the procedure until 2 + 2 k e > A gives the result. o 

Remark 4. As an example where the limit form is non-local, let us consider the case 
where (E^J 7 ) = (Q,A) and let 9 be a measurable map from f2 into itself preserving 
the probability P and defining a strongly mixing endomorphism A G A i— > 9~ 1 (A), i.e. 
such that 

HmE[fo9 n .g]=EfEg V/, 9 GL 2 (P) 

n 

then taking for Y the identity map and Y n = 9 n with a n — 1, hypotheses (HI) to (H3) 
hold onV = L°°(P) with ~A[tp\ = A[tp\ = A[tp\ = -tp + R[<p] and S[(p\ = varcp. o 

We introduce now the fourth bias operator^ defined under (HI) and (H2) on T> as 

)\= 1 -(A-A). 
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By E y ^[v?]x] = lim n E[(<p(Y n ) - ip(Y))( X (Y) + X (Y n ))/2] we see that 4 represent the 
asymptotic error from the point of view of an exterior observator according the same 
weight to both the theoretical and the practical models and measuring the error al- 
gebraically on the same axis. Because of the properties of proved below, )\ will be 
called the singular bias operator. 

We shall say that an operator B from V into L 2 (P Y ) is a first order operator if it 
satisfies 

B[<px]=B[<p] X + <pB\x] Vy9,xe© 

Proposition 1. Under (HI) to (H3) 

a) the theoretical variance lim n a n E[(<p(Y n ) — (p(Y)) 2 ip(Y)} and the practical vari- 
ance lim„ a n E[(ip(Y n ) — ip{Y)) 2 ip{Y n )} exist and we have V</?, X , G V 

hm n a n E[( V (Y n ) - V {Y)){ X {Y n ) - X {Y))^{Y)\ = E Y [-A\^] X + 4&]<PX ~ Mv\x^\ 
hm n a n E[( V (Y n ) - tp(Y))( X (Y n ) - X (Y))^(Y n )} = E Y [-A[<pi/>] X + A[i/>]<p X - A[<p] X i/>] 

b) These two variances coincide if and only if ){ is a first order operator, and then 
are equal to E Y [r[if\ip]. 

Proof. The part a) comes directly from the definition of A and A. The difference 
between the two expressions in ip, X , tp is 

2E Y [{A[^]-m\v-m^)x] 

and vanishes iff is first order. o 

A sufficient condition for the equality of the theoretical and the practical variances 
is given by 

Proposition 2. Under (HI) to (H3) // there is a real number p > 1 s.t. 

lima n E[( V (Y n ) - <p(Y)) 2 \ii(Y n ) - ^j(Y)\ p ] =0 G V 

n 

then 4 is first order. 

Proof. Let be A G [0,2[. We have \a n E[(ip(Y n ) - v?(F)) 2 (^(F„) - i/j(Y))]\ 

< a n E[\ V {Y n ) - V {Y)\ x \^{Y n ) - V {Y)\ 2 - x \^{Y n ) - ^{Y)\} 
the case p — lis obtained taking A = 0. If p > 1 we go on with A > 

< (aJL[(<p(Y n ) - <p(Y)fl)W(aJL[(<p(Y n ) - <p(Y))*MY n ) - ^Y)^ 2 ^])^/ 2 
the result follows taking 2/(2 — A) = p. o 
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In particular under (HI) to (H3), if the locality condition of theorem 1 is fulfilled 
then ^ is a first order operator. 

Remark 5. In the frequent cases where 4 is a first order operator, A captures all the 
diffusive part of A and of A. 

Similarly, we can remark that for deterministic approximations, the operator A is 
often nought (see prop. 17 in fine). For example let us consider the ordinary differential 
equation 

x t = x + / f(x s )y s ds 
Jo 

approximated by the Euler scheme 

Xt =x + / f(x^ ns]/n )y s ds 
Jo 

even if we suppose x to be random, errors are of deterministic nature and as soon as 
/ is C 1 with at most linear growth and f Q y 2 s ds < +oo we have for ip,x £ C\ (bounded 
with bounded derivative) 

nE[((p(x?) - ip(x t ))x(x t )} -> E[u t ip'(x t )x(x t )} 
nE[(<p(x t ) - <p(x?)) X (x?)] - -E[u tV '(x t ) X (x t )} 
and E[((p(x?) - ^{x t f} -> 

where u t is given by u t = —\ / * f'(x s )f(x s )y 2 s e^ f ' i - Xa ^ yada ds. (cf. [17] theorem 1.1). 
Thus 

A[(p](x) = E[u t ip'(x t )\x t = x} = -A[(p](x) 
and we have A = A and A = 0. o 

Let us derive some consequences of the only (HI) hypothesis. Under (HI) we may 
consider the symmetric positive bilinear form 

e[<p, X ] = -EL%]x + <pA\x] - A[<px}} 
Proposition 3. Under (HI) the following conditions are equivalent 

1) (m) 

2) (H3) 

3) (e, V) satisfies the following sufficient closability condition 

Xn^V Xn -* in L 2 =>- e[<p, Xn] ^0 V</? G V. 

Proof. Since lim n a n E[{y{Y) - <p(Y n )) X (Y n )] = E Y [A\x]<p-A[<px]], (H2) is equivalent 
to 

V</? e T> x l— > ^y[A[x]^ — A[ipx}} is continuous on V in L 2 
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which is equivalent to 

V</? G T> x M e VPi x] continuous on V in L 2 , 

i.e. equivalent to this bilinear form be continuous at which is the condition of the 
statement. o 

Proposition 4. Under (HI), the conditions of the preceding proposition are fulfilled if 
G T> 

Xp EV, Xp^OinL 2 => hm^ma n E[( Xp (Y n )-Xp(Y))ip(Y)] = 0. 

p n 

Proof. The condition of the statement means x ~^ ^ , [^[x] < p] continuous, i.e. T> C 
D(A ) hence by remark 1 hypothesis (H2) holds. o 

Remark 6. If V</? G V the conditional expectation a„E[(/?(y n ) — (p(Y)\Y = y] con- 
verges weakly in L 2 (¥ Y ) then (HI) is satisfied, because the weak limit is necessarily an 
element of L 2 . 

Proposition 5. Under (HI), if the law of the pair (Y n , Y) is asymptotically symmetric 
in the following sense : 

lima n E[ V (Y n )^(Y) - <p(Y)if>(Y n )] =0 Vy?,^ G V 

n 

then the conditions of proposition 3 are fulfilled, A — A = A and)\ = 0. 
Proof. Taking -0 = 1 gives lim„ a ri E[(p(Y n ) X (Y n ) — (p(Y)x(Y)] = hence 
lima n E[(^(y) - ip{Y n ))x{Y n )] = lima< n E[(^(y n ) - <p(Y)) X (Y)] 

n n 

and (H2) holds with A = A. o 

We now come back to the situation where the only assumption (H3) is supposed. 

Theorem 2. Under (H3). If the form (£,B) (cf. theorem 1) is local, then the principle 
of asymptotic error calculus is valid on 

v = {F(f l7 ...j p ) : fiev, FeC^r.R)} 

i.e. hm n a n E[(F(/ 1 (y n ), . . . , f p (Y n )) - F(f 1 (Y), f p (Y)) 2 ] 

= EyEU Wi. ■ ■ ■ , ■ ■ ■ , /p)r[/i, /;]]■ 

Demonstration, a) Let us first give the argument in the case p — 1. When the 
form is local, by remark 3, lim n a n E[|/(F n ) — f(Y) \ k ] =0 V/6D for any integer 



k > 3. Let F E C 1 (M, M), writing the finite increments formula F(y) — F(x) 
(y — x) Jq F'(x + t(y — x))dt we have 



a n E[(Fof(Y n )-Fof(Y)f] = a n E[(f(Y n )~f(Y)) 2 ( / F'(f(Y)+t(f(Y n )-f(Y)))dt) 2 ]. 

Jo 

Let Pk be a polynomial uniformly close to F' on the closed ball P(0, \\f\\oo), considering 

(/) = a n E[(f(Y n ) ~ f(Y)nj; P k (f(Y) + t(f(Y n ) - f{Y)))dtf] 

= a n E[(f(Y n ) - /(F)) 2 ( (P k (f(Y)) 2 + Ef =0 , m= i ^m(f(Y)Y(f(Y m ) - /(y))™) 
= a„E[(/(Y) - /(F„)) 2 ((P,(/(y;)) 2 + Ef=o, m =i ^m(f(Yn)Y(f(Y) - f(Y n )) ' 

we observe that (/) has same limit as 

Mf(Y n ) - f(Y)) 2 {(P k (f(Y))) 2 + (P k (f(Y n ))) 2 ) /2] 



Oir. 



which converges to Ey [rf/jP^/)] by theorem 1. 

Now, \a n E[(Fof(Y n )-Fof(Y)f]-E Y [F' 2 of T[f}}\ 



< a n E 



(f(Y n )-f(Y)f 



(!o l F'{f{Y)+t{f{Y n )-f{Y))dt) 2 

-(iiPk(f(Y)+t(f(Y n )-f(y))dty 



+ 



a n E 



(f(Y n ) ~ f(Y)) 2 (fiPk{f(Y)+t(f(Y n ) - f(Y))dt\ 



Ey[P 2 ofT[f]] 



+ |E y [PMr[/]]-E y [F' 2 o/r[/]]|. 



Since sup n a n E[(f (Y n ) — f (Y)) 2 ] < +oo, the first and the last terms may be made small 
uniformly in n by a suitable choice of k, as the second term goes to zero when n f +oo, 
the proof in complete in this case. 

b) In the general case the finite increments formula writes 

v pi 

F(yi, y p )-F(x u ...,x p ) = Y^y—Si) / P/(yi, • • • , Vi-i, Xi+t(yi-Xi),x i+1 , x p )dt. 

i=i J o 

The local property of the form implies 

k 

(2) lima n E[J]|/ i (y n )-/i(^)|] = VA, . . . , / fe e P VA; > 3 

i=i 

by Holder inequality E| nLi X il < Yli=i( E l\ X i\ k ]) 1/k ■ Then the P roof proceeds simi- 
larly, approximating the derivatives F[ by polynomials P^ on the ball B(0, maxj ||/i||oo) 
of W using (2) and the property VV, x G ^ 

iima B E[(/ i (y B )-/ i (y))(/ i (y B )-/ i (y)) fo>(r„)x(y) + vM)) /2] = Mm, z,-]^] 
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which is consequence of theorem 1. 
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Let us end this part by a remark concerning the transportation of the four bias 
operators by image (cf. also [8]). 

Remark 7. Let Y n be an approximation of Y satisfying (HI) to (H3) on the same 
dense algebra V of bounded functions with constants. Let $ be a map from (E, T) to 
(G, Ti) such that the algebra T><$> — {u : u o $ e V} be dense in L 2 (P$ D y). Let us put 
Z n = $oF n andZ = $oy. Then Z n and Z satisfy (HI) to (H3) with and the 
same sequence a n : 

A*[u](z) = E Y [A[u o$]\$ = z] 
Az[u](z)=E Y \A[uo$]\$ = z] 

and similar relations for A<$> and A^. The Dirichlet form associated with A<$> is the 
image by $ of the Dirichlet form associated with A, and the properties of images of 
Dirichlet forms (cf. [9] chapter V) apply (square field operator, locality). If A is first 
order, A^ is first order. 

II. Examples. 

II. 0. Preliminary example. This is not stricly speaking an example but a part 
of Dirichlet forms theory itself. Let (£,B) be a Dirichlet form on the Hilbert space 
L 2 (E ,J- ', m) where m is a probability measure and let (Pt) be the strongly continuous 
contraction semi-group associated with (£,D). 

Let us suppose that the quasi-regularity assumption is fulfilled so that we may 
construct a Markov process Y t with P t as transition semi-group (cf. [23] chapter IV 
§3), and let us suppose also that the domain T>A of the generator (A, T>A) contains 
an algebra T> of bounded functions with constants dense in L 2 . Then for / e T>, the 
approximate forms 

£t[f] = \<f- Ptf, f >LHm)= ^fm[(f(Y ) ~ f{Y t )f] 



do converge (increasingly) when t I to £ [/] = 
hypothesis (H3) is fulfilled. The form (£,D) 
theorem 1. 

Here, as easily seen, we have 



- - < Af, f > (cf. [13], [9] or [5]). Hence 
is an extension of that one provided by 



A[f] = A[f] = A[f] = A[f] V/GD 

and the operator A vanishes. The above properties of Dirichlet forms hold either for 
local or non-local forms. Since ^ | +oo we see that the hypothesis (H3) may be satis- 
fied with a n t +oo the limit form being nevertheless non-local (cf. e) of theorem 1). 
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Note. In the whole article the positive symmetric bilinear forms of real functions £[f, g] 
are extended to complex functions, not as Hermitian forms, but as bilinear forms with 
the same symbol S, in other words 

+ if 2, 9i + m] = /2] + i£[fi, 92} + iS[f 2 , 9i] - S\f 2 , 92}- 

II. 1. Error in the Glivenko-Cantelli theorem. Let us begin with a simple 
one-dimensional example related to the Glivenko-Cantelli theorem. Let X be a real 
random variable with continuous distribution function F anf let X k be i.i.d. copies of 
X. 



If we put 



n 



fc=l 



and Y = F(X), we have Y n — > Y a.s. and we may study Y n as approximation of Y. 
Thinking for instance X and the X^s to be simulated by the inversion method shows 
that the pair (Y n ,Y) has the same law as (- Y^=i ^-u k <u, U) where U is uniformly 
distributed on [0, 1] and the Uk are i.i.d. copies of U. Since our framework involves 
only the joint law of Y n and Y we may work with (Uk,U) instead of (Xk, X). Choosing 
a n = n and V = £{x 1— > e 2l7Tpx , p e Z}, we obtain that hypotheses (HI) to (H3) are 
fulfilled with 

%](y) = ^/(y) 

4M(v) = ^^'(v) + (i- W(v) 

m(y) = (y-W(y) 

the theoretical and practical variances coincide and D = {/ e £ 2 [0, 1] : /' in distribution 
sense e L^QO, 1[) and y 1— >• ^/y — y 2 f'(y) £ L 2 (dy)}. Let us give some indication 
on the proof of hypothesis (HI) for instance. Denoting the conditional law given 
Y = y, we have to study 



n 



f E y [(e 2iwpYn - e 2i7rpy )e 2i7vqy ]dy = n I ((e 2iwp/n y + 1 - y) n - e 2 ™ py )e 2i7Tqy dy 
Jo Jo 



which may be expanded as 

= n ^exp{-^-^(27rp) 2 }(l + e(n,y)) - l^j e 2l7T ^ y dy 

where e(n, y) goes to zero uniformly in y when n — > 00. Using e" A — 1 = —A e~ tx dt, 
the dominated convergence theorem applies and the limit is 

f - y -^(2np) 2 e 2 ^ y dy = E Y [^-^ V " (Y) X (Y)] 
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for ip(y) = exp2inpy and x{v) = exp2i7rgy. 

II. 2. Typical formulae of finite dimensional error calculus. 

II. 2. a. Let us consider a triplet of real random variables (Y, Z, T) and a real random 
variable G independent of (Y, Z, T) centered with variance one. We are interested in 
the approximation Y £ of Y given by 

(3) Y £ = Y + eZ + ^JeTG. 

In the multidimensional case, Y is with values in MP as Z, T is a p x g-matrix and G 
is independent of (Y,Z,T) with values in M, q , centered, square integrable, such that 
E[GiGj] = Sij. 

Operator A. 

Proposition 6. If Z and T are square integrable, if is C 2 bounded with bounded 
derivatives of first and second orders (ip e C 2 ) and if x is bounded, 

-E[(ip(Y £ ) - <p{Y)) X {Y)\ -> E Y [A[cp} x ] 

where A[<p](y) = E[Z\Y = y]<p'(y) + \E[T 2 \Y = y] V " \y) . 
In the multidimensional case 

A[ip\{y) = R[Z t \Y = y]Vtp(y) + ± £ E[(TT%\Y ^yj^y). 

ij 

Proof. Let us give the argument with the notation of the case q — p — 1. The 
Taylor-Lagrange formula applied up to second order gives 

±E[(^(Y £ ) - y?(Y)) x (Y)] = ELV(Y) X (Y)] 

+ \E[(eZ 2 + 2^/lZTG + T 2 G 2 ) £ ip"(Y + ab(eZ + ^TG))2adadb X (Y)\ 

(note that ZTG and T 2 G 2 G L 1 because of the independence) and this converges by 
dominated Lebesgue theorem to E[Z(p'(Y)x(Y)] + \E[T 2 V "(Y)x(Y)}. o 

Quadratic form and operator A. 

Proposition 7. If Z and T are square integrable, if if and x are C\ 

-E[( V (Y £ ) - <p(Y)(x(Y e ) - x(Y)} - E[TV'(Y)x'(Y)] 
and in the multidimensional case 

\e[{ V {Y £ ) - <p(Y)( X (Y £ ) - X (Y)\ -> E[(V^)*(Y)TTV X (Y)]. 
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Proof. The demonstration is similar with a first order expansion. o 

In order to exhibit the operator A, we must examine the conditions of an integra- 
tion by parts in the preceding limit. Let us put 9ij(y) = E[(TT*)jj|y = y] so that 

E[(v^)'(y)Tr*vx(r)] = E« E y[^«x / i]- 

Proposition 8. If Z and T are square integrable, if fori,] = l,...,p the measure 
OijWy on MP possesses a partial derivative in the sense of distributions 9j(%Py) which 
is a bounded measure absolutely continuous w.r. to Py 7 say p^-Py, then as soon as 
9ij and pij G L 2 (Py) the form £[<p,x] = \ Ey[</^%Xj] ^ s closable on the algebra 
V = Cl, hypotheses (HI) to (H3) are fulfilled and 



ij ij 



Proof. We have 



ij ij 

and the equality 

J OijdMxWY = ~ J <p'iXPijdW Y 

valid for ip, x £ extends, under the assumptions of the statement, to ip, x £ C\. 
This yields 

\ WAjX' 3 ] = -\J(E 9 ^ + E Pa^xdPy. 

ij ij ij 

Q.E.D.o 

The operator A depends only on T, not on Z. We obtain A by difference : 

m<p\ = \J2 9i w"j + EE - z M 

ij j i 

where Zj(y) —~E[Zj\Y — y\. At last, A is first order : 

= J2(zj ~\^PijWj- 

j i 

Remark 9. The results of this section II. 2. a) would be identical with an approximation 
of the form 

(4) Y £ = Y + eZ + T.B £ 

where B is a centered Brownian motion vanishing at zero independent of (Y, Z, T) since 
only the joint law of (Y, Y £ ) is used. 
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The question is very close to the classical approach of Kolmogorov [20] to study- 
Markov processes starting from the assumptions 

lim^o M( X t+h - X t )\H = b(t,X t ) 

\im hl0 \E[{X t+h - X t ) 2 \F t ] = a(t,X t ) 

Indeed, it is easy to see that the representation (4) occurs naturally for Ito processes 
and for diffusion processes given by an Ito equation : 

On a filtered probability space (jF t ), P), let B t be an (jF t )-Brownian motion 
centered vanishing at zero and let £ be an Ito process defined by 

it = £o + / o- s dB s + b s ds 
Jo Jo 

where the processes a and b are adapted and continuous at zero in L 2 (P) and L X (P) 
respectively. Then approximating £ by £ 4 is equivalent to approximating £ by £ t = 
£o + tb + °~oBt because for y? G C 2 

lim t _ >0 iE[( V 9(6)-^(6) 2 |^o]=0 

lim^o \m<p(St)-<p(l)\Fo]=0 

as soon as for instance Eof is bounded in a neighborhood of zero, as seen by application 
of Ito formula and standard inequalities. 

II. 2. b. Series with independent increments. 

Let be 

s = + ^ 

n=l 

where X n , Z n e L 2+e , Z n centered, (X n , Z n ) i.i.d., we approximate S by its partial sum 

C _ V™ X k I Z k 
° n ~ l^k=l k 2 " r k ' 

Using Burkholder inequality, we observe that nE[|S f — S , „| 2+£ ] — > as n — > oo. Thus, 
taking V = C|?, we have for ip, x £ ^ 

limnE^S) - y^)) 2 ] = limnE[(S - S n )V 2 (S n )] = E[Z 2 ]E[y/ 2 (S)] 

n n 

\imnE[(<p(S) - <p(S n )) X (S n )} = limnE[(S - S n )<p' (S n ) X (S n ) + hs - S n ) 2 <p" (S n ) X (S n )} 

n n Z 

= ^E[Z 2 W(S) X (S)] + E[X 1 ]E[<p'(S) X (S)}. 
We can conclude that hypothesis (H2) is satisfied and 

am = + nxiW- 
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Assumption (H3) holds as soon as the law of S satisfies the Hamza condition ([13] 
p. 105) and then the Dirichlet form is local. 



II. 2. c. Tails of martingales. 

Let us first consider the classical case of Polya's urn in its simplest configuration 
with two colors, one ball added each time, and an initial composition of one white ball 
and one black ball. 

The ratio X n of white balls after the n-th drawing satisfies 

X n+1 (n + 3) = X n (n + 1) + l Un+1 <x n 

where U n+i is a random variable uniformly distributed on [0, 1] independent of T n = 
a(X , . . . ,X n ), i.e. 

X n+ i = X n H —(lu n+1 <x n - X n ). 

n + 3 

Let X be the (a.s. and L p , 1 < p < +oo) limit of the bounded martingale X n , we 
study the approximation of X by X n . 

We note that lim n nE[|X — X n | 3 ] = as easily seen using Burkholder inequality. 
Then, taking for T> the functions of class C 3 on [0, 1] vanishing at and 1, we have 

\imnE[(<p(X)-<p(X n )) X (X n )} = \imnE[{X-X n ) V \X n ) X {X n ) + hx-X n ) 2 ^' {X n ) X {X n )} 

n n Z 

and 

E(X-X n )} -E^ y—y 2 

k=n 

6n 

because X n X and X is uniformly distributed on [0, 1], as easily verified. 
We obtain 

limnE[(y?(X) - tp(X n )) X (X n )} = ±-E[<p" (X) X (X)] 

n IZ 

limnE[(^(X) - <p(X n ) 2 ] = l -E[p' 2 {X)\. 

n <Q 

Hence (HI) to (H3) are fulfilled A[<p] = and A[(p] = ±(p" so that A = A and 

A = 0. The limit error structure is the uniform error structure on [0, 1]. This analysis 
could be easily extended to any configuration of Polya's urn, mutatis mutandis. 

More generally, this kind of asymptotic behavior appears, under regularity assump- 
tions, for the approximation between a martingale and its limit. 
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Let M n = Y^i=iXi be a martingale w.r. to the filtration T n . Let us suppose M n 
centered, square integrable s.t. YliLi ^-^i < +°° an d ^ us P ut a n = Y^L n +i ^-^i- 



Proposition 9. Supposing ^- sup i>ra \Xj\ — > m probability, E[^sup i>ra X?] bounded 

in n, E[(^(^ i>n Xj) 2 ) p ] bounded in n for some p > 1, and ^ X]j>n^? — * C 2 
probability, then 

a) y- J2i> n Xi ==> where Z has for characteristic function E[e~^ 2 ' 2 ] 
^ hypothesis (H2) is satisfied andWip,x £ 

lim n ^E[( V 9(M 00 ) - <^(M n ))x(M n )] = |EZ 2 E[ V 9 // (M 00 ) X (M 00 )] w/mi awes 

AW\ = \nz 2 ] <p". 

c) lim n 4E[(^(Moo) - (f(M n )) 2 ] = EfZ^Er^fMoo)], hypothesis (H3) zs satisfied 
as soon as the law of satisfies the Hamza condition, then A — A = A and A = 0. 

Proof. The proposition is a direct consequence of a result of Hall and Heyde ([14] §3.5 
p. 76 et seq.). 

II. 3. Conditionally Gaussian case 

Let us begin with the finite dimensional case before applying the approach to pro- 
cesses. 

II. 3. a. Let Y be a r.v. with values in V = (V^)^^...^ be an application 
from M. d into symmetric positive d x d-matrices, £j be r.v. with values in M. d which 
conditionally given Y = y are i.i.d. Gaussian with common law J\f d (y,V(y)). We 
consider 



as approximation of Y. We take V = C{x i— > e i< "' a:> , w G M d } and a n = n. 



Lemma 1. Let us suppose E[trace(V(Y))] < +oo 7 then 

1) lim n nE[(^(y n ) - V {Y)f] = EyEy=i </^d Vy? G X>, 

,2) £/ie hypothesis (HI) zs satisfied and A[ip\(y) = | X^ J=1 V i j(y) l Pij(y) V<£> G £>. 

The proof proceeds without difficulties as in the preceding example by finite expan- 
sions of the exponential functions. 

Lemma 2. Le£ suppose E[trace(l / (F))] < +oo and £/ie following regularity con- 
dition: the measure V^Fy has a partial derivative ^(V^Py) in the sense of 
distributions which is a bounded measure absolutely continuous w.r. to F Y , say Pijfy, 
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then as soon as Vij and pij G L 2 (Fy), (HI) to (H3) are fulfilled and 

= Part 

Proof. The condition of the statement allows to perform an integration by parts in 
the limit obtained in lemma 2. That gives (H3) hence (H2) as well. o 

Remark 9. There are several sufficient conditions in order that the form £[u, v] = 
EyEij«H^-] be closable 011 C%(R d ) (cf. [13] chapter 3 §3.1 and [23] chapter II §2). 
Suppose such a condition holds, then by the argument of remark 2 the hypothesis 
(H3) is equivalent to T> C VA where (A, VA) is the generator of the smallest closed 
extension of (S,C%(R d )). 

II. 3. b. Let us apply this to the approximation of processes. Let Y be a real 
process indexed by a set T. Let us consider Y as a measurable map from (Q, A, P) 
into (E,F) = (R T , (B(R))® T ). Let y t be the coordinate mappings from E on R. We 
consider the algebra 

V = £{e iuiyt i + - +iUkyt K, Uj G R,tj G T}. 

Thanks to the monotone class theorem, T> is a dense algebra in L 2 (E, J 7 , Py). We put 
Y t = y t oY. 

Let = (Ci)teT be a sequence of real processes such that, conditionally given Y 
the £ J are independent with the same Gaussian law with 

ml\Y]=Y t 

E[(£-Y.)(£-Yt)\Y]=C., t (Y a ,Y t ) 

the function C Sjt (x 1 , x 2 ) and the process Y being such that E[C Sjt (Y s , F t )] < +oo Vs, t G 
T. 

We approximate F by the process Y n = \Y^j=i& ■ The results obtained in the 
finite dimensional case give the following proposition : 

Proposition 10. If the marginal laws of Y P(y tl , ...,Y t )(dyt 1 ■ ■ -dyt k ) possess partial 
derivatives in the sense of ditributions which are bounded measures absolutely con- 
tinuous w.r. to P(y tlr ..,y t ), say \j{y tl , . . . , y tk )Py ; then hypotheses (HI) to (H3) are 
verified and for <p eD we have 

A[<p](y tl ,...,yt k ) = \ E?j=i c u,ti (vu , Vtj ) dy ^o ytj 

• • • >z/t fc ) — I Eij=i Ct i ,t j (yt i ,yt j ) d yt .d yt . + Eij \j(yti, ■ ■ ■ >ytj^f: 
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Under the hypotheses of proposition 10, the form (£, D) is local and theorem 2 on 
asymptotic error calculus applies. Let us also remark that if we define the operator # 
on V by 

j=i yt J 

when ip depends only on y tl , . . . , y tk , we have 

E[(<p#) 2 \Y] = T[<p](Y) 

and f ((/9 # ) 2 rfPy(y) = 2£[ip]. It follows that the operator # extends uniquely to D in a 
closed operator satisfying for F e C 1 fl Lip(IR p ,IR) 

(F(f 1: ...J p ))# = J2F;(f 1 ,...J p )f* V/i,...,/ p GD 

In other words # plays the role of a gradient w.r. to the Dirichlet form (£,D). 

Special case 1. F is a real process indexed by R + , and £| = F t + W/ where W/ 
are independent standard Brownian motions independent of Y. If the marginal laws 
of Y have densities f tl ,...,t k s.t. dft g^" tk = \{yt^ ■ ■ ■ ,Vt k )ft u ...,t k with A* G L 1 (R fc ), the 
construction applies and yields a Dirichlet form with square field operator Y s.t. 

r[<p](Y) = -^UiUjU M,ip 2 (Y) for <p(Y) = e^" fc5 V 

ij 

Suppose Y possesses second order moments, then the linear forms £ = Ylp=i a p0^t p+1 — 
Y tp ) are in D and F[£] = Ylp=i fl p(^+i — t p ), so that on step functions / 

T[J fdy] = J f(s)ds. 

Thus this error structure may be called the Ornstein- Uhlenbeck structure on the process 
Y. 

Special case 2. With the above notation, let us assume that 

a = Y t - h(t)wi 

where h is a deterministic function. With the same hypotheses as above, we have on 
step functions / 

T[J fdy] = n(J f(t)d(h(t)W t )) 2 ] 
Special case 3. Suppose eventually 

il =Y+ f h(s) dWi 
Jo 
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then Ct u tj = Jq M3 h 2 (s) ds and T[J f dy] = J f 2 (s)h 2 (s) ds. We obtain a structure 
which may be called the Ornstein-Uhlenbeck structure with weight h on the process Y . 

II. 4. Natural inaccuracy of the Brownian motion simulated by the Donsker 
theorem. 

We begin with the simplest case of one dimensional marginal laws which is here 
nothing else than the central limit theorem. 

II. 4. a. Natural inaccuracy in the central limit theorem. 

Let be S p = YTi=i ^ where the randon variables Vi are i.i.d. centered with variance 
a 2 . We consider two indices m and n linked by the relation 

n = n{m) = m + k{m) with 9y/m < k(m) < -y/m 

9 

for a 9 e]0, 1[. 

Let us observe some evident properties : lmv^oo ^ = 1; G [0, |]; y/n > 

^VTT9>9VTT9(n-m);^-^< 1 t. 

Writing -^=S n - ^S m = -^=(S n - S m ) + - ^=)S m and using these properties, 
shows that A=S n 7=S m — > a.s. 

We consider the mutual approximation of -^S m and -j^S n (which is an obvious 
extension of the framework of part I). That is for A we study 

(5) a(m)E[M-Ls m ) - ^( J=5 B ))x(4=5 B )] 
and for A we study 

(6) a(ro)E[M-^S B ) - ip{^=S m )) X {-^S m )] 



with ct(m) = (so that 9y/m < a(m) < ^y/rn). For the algebra V we take the 
linear combinations of imaginary exponentials. 

Proposition 11. Suppose the Vi 's possess a third order moment, then hypotheses (HI) 
to (H3) are fulfilled and for <p &T> 

A[<p](x) = A[<p](x) = A[<p](x) = \o\" - \x^'. 

The Dirichlet form is the Ornstein- Uhlenbeck form onM, ( endowed with the normal law 
M^a 2 )). 
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Proof. For A, taking (p(x) = e lux , x{ x ) = e%m we have to look at 

J m = a(m)E[e i( ^ + i^ )Sm+ ^ {Sn ~ Sm) - e'W 5 "]. 
Let £(t) = E[e** yi ] be the characteristic function of the V^'s 

Jm . a(m)K( _|. + ^ m ^ )r ™_ K( ^ )r] 

using the existence of a third moment we have 

(7) ]ogt(t) = ~o*(l + at + to(l)) 

CT 2 CT 2 

and this allows to write J m = a(m)[e~^~ Km — e~^~ Lm ] with 



& a /TO/ ' 

K m = (u + vf + -t=(u + w) 3 + 2w(-^= - 1) + -= o(l) 



and L m = (u + v) 2 + -2=( u + vf + ^o(l). This yields 

J m = aHe-^K-^M^-l) + ^o(l))] 

hence 

hm J m = i^W-^W = ~ / MdAffro*) 

m 2 2 J R 

what gives easily the proposition. o 

II. 4. b. The Donsker case. 

Let the V?s be as before and 



, / [«*] 
X n (i) = — VV fc + (ni 



[n*]Vj nt ]+i 



for t G [0, 1], [nt] denoting the entire part of nt. 

The laws of the variables X n are probability measures on C([0, 1]) as well as their 
limit in law which is a Brownian motion by Donsker theorem. 

For the algebra V we take the linear combinations of exponential of the form ip(X) = 
exp {iX(f)} where X(f) = f(s) dX(s) and with / e C 1 in order that £ /(s) rfX(s) 

may be defined as X(l)f(l) — X(0)f(0) — J* X(s)df(s) for the general coordinate 
process X(s) on C[0, 1]. As easily seen the algebra V is dense in L 2 (C([0, 1]), //) /i being 
the Wiener measure. 



21 



Thus we have X n (f) = v^X]fc=o /fc/r^ 1 ^™ f( s ) ^ s Vk+i an d for studying the operator 



A we have to look at 



M m = a(m)E[( V (X m ) - <p(X n )) X (X n )] = a(m)E[(e lX ^ - e iX ^)e iX ^\ 



aimW, 



fc + 1 fc | 1 \ fc | 1 

i>/^Zk£j k m f(s)dsV k+1 iV^EkZoh" f(s)dsV k+1 \ iV^EkZoh 71 9(s)dsV k+1 



We take as before a(m) = m/k(m). 

Proposition 12. Suppose the V* 's possess a third order moment, then hypothe- 
ses (HI) to (H3) are fulfilled. We have A = A = A on T>. The Dirichlet form 
is the Ornstein-Uhlenbeck form on the Wiener space (with a Brownian motion s.t. 
< B > t = a 2 t) normalized so that the square field operator satisfies T[J^ h(s) dB s ] = 
^h 2 {s)a 2 ds V/i E L 2 ([0, 1]). 

Since the Dirichlet form is local, some limits are automatically obtained (theorem 
2). Since A = 0, the theoretical and practical variances coincide (prop. 1). 
Proof. For studying A we consider the quantity M m defined above. By the third 
moment assumption, the characteristic function £(£) of the V$'s satisfies (7) and we can 
write 



M m = a(m) 

(8) 
with 



n— 1 



\_)=m 



2+1 m-l 



n-l 



j=m 



3 + 1 



n k=0 

= a(m) 



e 2 



e 2 



n-l 



j=m 



3 + 1 



3 + 1 



N m = Y,<] 3 g) 2 + anyfti(j. g) 3 (l + o(l)) 

n n 

m-l r fc+1 r h±l fe+i fc+i 

+ ^(v^ j k m f + V^ j k " g) 2 + a(V^ f + V^ J k " g?d + o{i)) 

k = m n m n 

En „ In „ 

n( y. (f + g)) 2 + anV^(J 3 (f + g)) 3 (l+o(l)). 

k = n n 

J + 1 J+1 

Using mY^Joili^ ff = Jo f 2 + and su Pi Iv^jT^ #1 < ^IblU we obtain 

1 1 m l fc+1 fc+1 ^ 

iV m = / s 2 +/ f 2 + 2j2V^V^ [ m f [ " g+^ [ (f + gf(l + o(l)) + -0(l) 
Jo Jo J± J± V™Jo n 



n-l 
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and 

fi 



m 



(/ + </)' + -7= / (f + 9Y(1 + o(l)) + -0(1) 
o v m Jo n 



Putting these expressions in (8) leads to 

M m = « W («p[-£ jf (/ + ,)«]) [(-^K^f - 1) jf ' /, + ^o(l) 
Eventually, for y?(X) = exp[i Jjj 1 /dX] and x(X) = exp[i gdX] we get 

i)E[(<p(X m ) - <p(X n )) X {X n )\ - y exp[-y j\f + gf ds] jf ' fg ds 



aim) 



In order to recognize the obtained limit, let r ou be the Ornstein-Uhlenbeck square field 
operator on the standard Wiener space (s.t. T ou [J hdB] = J h 2 ds). We have by the 
functional calculus 

Er OT [e^o^, e l /o^] = - j fgdsE[e^^ dB ] = - j fgdsexp[-^ J (/ + gf ds}. 

It follows that for a Wiener measure s.t. < B > t = u 2 t and an Ornstein-Uhlenbeck 
structure s.t. T[Jq hdB] = J h 2 d < B > whose generator will be denoted (A, VA), we 
have V C VA and 

<Aif, x > =-£[<p,x] = -m r [^x]] 

= t/o 1 f9dsexp[-^-J^(f + g) 2 ds] =< Atp, x > ■ 

The operator A is therefore symmetric on T>, which implies A = A — A and the propo- 
sition is proved. o 

Comment. As noted already by Louis Bachelier, assets quoted on the organized 
markets look like Brownian paths. This is displayed in any course in mathematical 
finance in order to introduce modelling by diffusion processes and stochastic calcu- 
lus. Nevertheless some concrete facts prevent this ressemblance from being accurate at 
microscopic scale. First because the spot is only defined at discrete instants, second be- 
cause a Brownian path possesses ideal properties (like the fact that it cuts uncountably 
many times every level that it reaches) that cannot be verified by material recordings. 

So that in order to be completely pragmatic, we might replace in financial models 
any Brownian motion by a random walk by application of Donsker theorem with n 
sufficiently large and consider the Brownian motion of the model is nothing else than a 
class of such sufficiently fine random walks. An infinite precision for stochastic calculus 
in finance is therefore a priori absurd and, by the results of this section, we may repre- 
sent the intrinsic fuzzyness of these computations by the Ornstein-Uhlenbeck form on 
the Wiener space. This is a justification of the approach proposed in [4]. 
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II. 5. Empirical laws and natural inaccuracy of the Brownian bridge 



If (V n ) are i.i.d. real random variables < V n < 1 with distribution function F and if 
F n (x) = - Y17=i ^-Vi<x is the empirical distribution function, then ^/n(F n — F) converges 
in law on the Skorohod space to a transformed Brownian bridge Bp( x ) ~ F(x)B 1 (see 
for instance [3] p. 141). 

Considering the V^'s are simulated by the inversion method shows that this result is 
a consequence of the special case where the V^'s are uniformly ditributed on [0, 1]. From 
now on, we restrict to this case. Putting Z n (x) = y/n(F n (x) — x) we are interested in 
the limit R m = a(m)'E[(ip(Z m ) — ip(Z n ))x{Z n )} for m and n linked as in the preceding 
example, with also a(m) = m/k{m) and for if, x G V where V is the algebra of linear 
combinations of imaginary exponentials of the form ip{Z) = exp{i J* f(s)dZ(s)} = 
exp{-ifi Z{s)df{s)} for / G ^([0, 1]). 

Thus 



Rn 



a(m)E 



=i(/(Vk)-//) _ e ^Efc=i(/(v fc )-//) 



Putting f = f — Jq fds and g = g — gds and denoting 77 and ( the characteristic 
functions of g and f + g gives 



aim) 



(p(m)) m ( (v(^)) n - m 



- (C(— ))' 

'n 



with p(m) = Wyj ^~ 9 \. 
The estimates 



P {m) = 1 - ^e[(/ + + (^nn(f + + 

C(^) = 1 - £E[(/ + ~g?] + (^) 3 |E[(/ + ~gf\ + 4,0(1) 

v(^) = 1 - m 2 \ + (^) W] + &o(i) 



allow to obtain 



\im m R m = e -^[(/+^) 2 ](-iE[5 2 ] + §aw(<7, / + </)) 

= I e -| E [(/+5) 2 ] COt ;(/ )5 ). 

In order to recognize the limit, let T ou be as before the Ornstein-Uhlenbeck square field 
operator on the Wiener space and Z t — B t — tB 1 . We have 

T ou [e l ^ dZ ,e l ^ dZ ] = -e^^ dZ T ou [ J fdZ, J gdZ] = - e *J(/+^( J fgds-J fds J gds) 



and consequently 



l e -hif(f+9) 2 ds-(f(f+g)ds)Z] 



lim m R 

n 



(J fgds - j fds j gds) 
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As before denoting A the Ornstein-Uhlenbeck operator, we see that A = A on V and 
therefore A = A — A, the Dirichlet form is the image of the Ornstein-Uhlenbeck form 
on the Wiener space. 

II. 6. Erroneous empirical laws and generalized Mehler type structures on 
the Brownian bridge 

We still consider a sequence V = (V)j 6 N of i.i.d. random variables uniformly dis- 
tributed on [0, 1] and the empirical distribution function 



I n 

f„(x) = ~y^iv<x 

II < * 



n . 



but the problem that we tackle is different. We suppose that the law of V 1 is not 
perfectly known. We assume that there is a sequence of r.v. approximating V 1 
and copies (U % m) V 1 ) of (U^V 1 ) such that the sequence (U l m) V % )i & \ be i.i.d. and we 
suppose that the law of has support in [0, 1] with distribution function F m . We 
define the emprirical distribution function 



We are interested in the approximation 

1 n 

of the process 



n . 



1 

z n = —= y (iv*<x — x). 

J=l 

We take V = C{^{z) = exp[-i J* z(t)df(t)}, f G C\% 1])} and we study 

Q m , n = PmVMZ?) - V (Z n )) 2 } = /3 m E[(exp[i— £/(C/4)] - exp[*— £ /(V% 2 ] 

v 1=1 v 1=1 

in which the symbol ~ represents the centering operation. 

Denoting rj [resp. rj m ] the characteristic function of f{V x ) [resp. f{U l m )}, 9 m the 
characteristic function of f(U^) + /(V 1 ), a 2 [resp. cr^J the variance of / [resp. of 

f( u L)h we have 

Qm,n = (3 m [(Vm(^)T ~ 2(fl m (^=)) B + M^=))1 
'71 \/Tl \/Tl 
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The estimates 



ri m (t) = l-^+^o(l) 

0m(t) = 1 - §E[(/(Eft) + 7(V%] + t 2 o(l) 

^(i) = l-^(T 2 +t 2 0(l) 



give 



Qm,n /3m& 



- 2cr 



e -2( CTm - CT 2 )+o(l) _ 2e 2 ( r 2 -lE[(/(^)+/(V 1 )) 2 ]+o(l) + expo (l) 



Remarking thatj^ 2 -±E[(/(£ft) +/(^)) 2 ] = §E[(/(Eft) + /(V 1 )) 2 ] + a 2 - a 2 m , if we 
assume E[(/([/^J — /(V 1 )) 2 ] ^ m -*oo and a 2 — o 2 m — > 0, we obtain 

Q™,n = /3 m e- 2CT2 (-E[(/(^J - fiV^)) 2 ] + o(l)) 



we can state 



Proposition 13. // there is a sequence (3 m — > +oo s.£. 

AnE[(/(^) - /(n) 2 ] - ^[/] 

where e[] is a quadratic form defined onC 1 ([0, 1]) closable in L 2 ([0, 1]) (with a Dirichlet 

extension non necessarily local, cf. remark 2) and supposing K[(f (U^)) 2 } — > E[(/(l /1 )) 2 ] 
then 

lim ArJEK^J" ^ - e^") 2 ] = V^ef/], 

hypothesis (H3) fulfilled and the limit Dirichlet form is the image by the Brownian 
bridge of the generalized Mehler type form on the Wiener space associated with the 
form e[.] (cf. [5] chapter VI §2.5 pll3 et seq). 

Proof. The hypotheses of the proposition imply what we needed during the above 
computation. It suffices therefore, as before, to recognize the limit as a closable form. 
But that comes from the functional calculus and the fact that the generalized Mehler 
type structure associated with the form e[\ satisfies T[J^ f(s) dB s ] = e[f\. o 
For example if — ^ Y^=i^Y l + a(V 1 )G p ) where the G v are i.i.d. reduced normal 
variables independent of V 1 and where a is continuous, then by the lemma 1 of the 
conditionally Gaussian case and the Hamza condition, the form 

e[f\ = IimnE[(/(tS - f(V 1 )) 2 } = f f'\x) a 2 (x) dx 

Jo 

is closable and the proposition applies. This generalized Mehler type structure satisfies 

(9) Y[f hdB]= f f' 2 {s)ds \/feH\[0,l]) 

Jo Jo 
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this structure may be constructed elementarily as in [5] or by the white noise theory. 



Remark 10. There exists an extension of Donsker theorem to the case where the 
variables (notation of section II. 4. b) are erroneous with a functional convergence 
in the sense of Dirichlet forms (see [6] and [10]). The limit structure obtained is the 
Ornstein-Uhlenbeck structure. This is related to the fact that the perturbation involved 
in this approach is a transversal perturbation of the random walk hence at the limit a 
transversal perturbation of the Brownian path (we will display this result in terms of 
an approximation procedure in section II. 7 below). 

Here instead, we change the law of the starting random variables : In the expression 
ofZ™ 

n 

if we simulate by the inversion method we see that has same law as 
-= 5><f»)-i<*)<. - F m (x)) = -= X)(1*<f»(.) - 

3=1 3=1 

where (X- 7 )^^ is a copy of (V^j^. We see that when F m changes, the path is longi- 
tudinally perturbed and so is the limit Brownian bridge. This explains a formula like 
(9). 

II. 7. Erroneous random walk and Donsker theorem. 

This example displays many similarities with examples II. 4, II. 5 and II. 6. We give 
only the framework and the results. 

Let U 1 be a centered square integrable r.v. approximated by also centered and 
square integrable. We suppose 

a m E[(U m - Uf] -> A. 
Considering i.i.d. copies (U % m) U l ) of (U^U 1 ) we look at 

[nt] 

* i=l 
[nt] 

X n (t) = -j=C£, Ul + ( nt "" [nt])U^ +1 ) 
v i=i 

and we study T mi „ = a m E[((p(X r ^) - (f(X n )) 2 ] for tp belonging to 

V = C{exp[i[ f(s)dX s ], f G C°° with support in ]0, 1[}. 
Jo 
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Putting a 2 = var(U l ) we find that 

lim T m>B = -A(/ fds)e- 2 ° 2 $f 2ds 

The limit Dirichlet form is once more the Ornstein-Uhlenbeck form on the Wiener 
space s.t. < B > t = a 2 t and T^ 1 / dB] = A f 2 (s) d < B > s . 

Comment. The Dirichlet- version of the Donsker theorem proved in [6] supposes the 
r.v. U 1 has a regular law allowing to carry a (non zero) Dirichlet form. This excludes 
the case of a discrete law. Here instead, we do not need such a restriction and the 
present construction applies for instance to the Cox-Ross-Rubinstein model approxi- 
mating the Black-Scholes model. But the convergence here is weaker than the one used 
in [6]. 

II. 8. Approximation of the Brownian motion defined through the Wiener 
integral, i.e. as centered orthogonal measure. 

Let X\ be a centered square integrable real r.v. approximating the variable X 1 
which is reduced Gaussian. Let (X^,X l ) ieN * be i.i.d. copies of (X^X 1 ). 

We assume X\ and X 1 satisfy the hypothesis (H3) with a n and an algebra T> 
dense in L 2 (A/"(0, 1)) of bounded functions containing the constants and the imaginary 
exponentials: Vy?, x £ 

(10) a n EMK) ~ V>PO)(xPO - xiX 1 ))] - A / v'x'dM^ 1). 

Jul 

Let be an orthonormal basis of L 2 (Ei, T\, fii) where (E 1 , Hi) is a cr-finite mea- 
sured space and let us consider the mappings 

oo 

/ G L 2 (Ei, JF i, fii) h J2<f>^ >X n 

k=l 

oo 

k=l 

We consider J n as an approximation of J and for the algebra V we choose 

V = £{$ : G CORM ^ $(0) = e i<f ' e> 

where / G L 2 (Ei, T\, fii) has a finite expansion on the basis 

here CORM denotes the set of centered orthogonal random measures on L 2 (E\,T\, fii). 
We study the limit of A n = a n E[$(J n ) - $(J)) 2 ] for $(0) = e i<f ' 9> with / = 
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We have 

An = Oi n (n?=i Ee 2 ^ - 2 n?=i Ee^(^+ Xl ) + n?=i Ee 2 ^* 1 

= «n (ll?=i Ee 2 ^ x « - n?=i Ee^( x - +Xl ) + n?=i Ee 2 ^ xl - f[? = i Ee^( x -+ x1 )) 
= a„(A n + B n ) 

and we may write 

Q q-l Q 

A n = l[^e 2lfkX " (Ee 2lf « x - - Ee 2 ^ x - +x1 ^ f[ Ee 2 ^ x - +X ^ 

q=l k=l k=q+l 

because the intermediate terms cancel and remain only the first and the last ones. 

Let us assume now in addition that the pair (X*, X 1 ) converges in law to (X 1 , X 1 ). 
Then 

\im n a n (A n + B n ) = £ J] Ee 2 ^ 1 limE[(e^ x « - e if " x ') 2 \ f[ Ee 2tf " xl 

g=l k=l k=q+l 

- EKe^-e^ 1 ) 2 ] 



= yE[e 2i ^=i^]liin 

' ^ n. 



a. 



E\e 2i ^ x 

q=l 

and by the assumption (10) this is nothing else than 

Q 

lim A n = E[e 2i ^f^ k ] V (-A/ 2 ) = E[e 2 * J ^ 

q=l 

We recognize once more the Ornstein-Uhlenbeck structure on the abstract Wiener space 
defined by J. We can state 

Proposition 14. If X* is an approximation of X 1 satisfying (10) and if (X^X 1 ) 4> 
(X^X 1 ) then the approximation J n of the centered orthogonal random measure J sat- 
isfies (H3) on V and yields the Ornstein- Uhlenbeck form. 

We would prove easily following the same lines that if the construction is done with 
different speeds for the different approximations of X k by X*, for example replacing 
X k n by X k defined by 

K = Kin) with -> a k , 

«4(») 

we would have 

a n E[( V (X k ) - V (X k )) 2 } ^ a k f ip' 2 dN '(0,1) 



and we would obtain the generalized Mehler type error structure on the abstract Wiener 
space defined by J associated with the quadratic form 



e[f] = J2 a of< 
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i.e. associated with the semi-group 

ptf = J2<f^« >e ~ aqt ^ 

(cf. [5] pll3 et seq). 

II. 9. Approximation of a Poisson point process. 

Let X be a r.v. with values in a metric space E endowed with its Borel a-field T . 

Let X n be an approximation of X satisfying hypothesis (H3) with the sequence a n 
and an algebra V Q composed of bounded continuous functions (containing the constants 
and dense in L 2 (F X ))- We suppose that the Dirichlet form defined by 

lima n E[(v?(X n )-v?(X)) 2 ] <p E V 

n 

is local. We denote by 7[.] its square field operator. 

Let fi be the law of X on (E, T\ Let (X^, X 7 ) be i.i.d. copies of (X n , X) and let J 
be an integer valued r.v. with Poisson law of parameter 1 independent of the sequence 

We consider the Poisson point processes 

j j 

3=1 j=l 

(Xa meaning zero). N n and N are r.v. with values in the space of point measures 
M. p on (E, T) equipped with the smallest a-field making all maps A i— > v(A) A e T 
measurable for v e M. p . 
We consider the algebra 

Lemma 3. V is dense in L 2 {¥ x ). 

Proof. By the chaos decomposition, it_is enough to prove that the constants and the 
elements of L 2 (F N ) of the form J fi dN ■ ■ ■ J f k dN where the functions f 1 ,...,f k are 
measurable bounded on (E,^) and where N denotes N — fi, may be approached by 
elements of T>. Since the constants are in V, it suffices to reach J fi dN ■ ■ ■ j fk dN. 
Now \im\^ (e x f ^ dAr — 1)/A = J ipdN and this gives easily the lemma. o 

We study the approximation of N by N n by looking at E n = a n E[($(iV n ) - $(iV)) 2 ] 
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for $(z/) = exp[i f ipdu}. We may write 

j j j 

E n = a n E[exp(2* ^ <p(X> n )) - 2 exp(* ^(Xj) + ^)) + ex P( 2 * E 



OO 1 

e 



= « n V" ^[(Ee 2l ^ x ^)P - 2 (Ee l ^ x ^ (xl ») p + (Ee 2 ^ xl )) p ] 

p=0 y 
oo _i 

p=o ^' 

OO _1 

p=0 

We apply a similar idea to what we have done in section II. 7 writing A p under the form 



A p = E 



n P 



fe=i 



Making now the additional assumption that (X^X 1 ) 4» (X^X 1 ) and using the fact 
that if G V Q is continuous and bounded, we see that a n (A p + B p ) has the same limit 
as 

aJEKe*^) - e^^XEe 2 ^ 1 ))^ 1 

but, since the form lim n a n E[(</?(X n ) — y?(X)) 2 ] is local, theorem 2 on asymptotic error 
calculus applies and gives 

lima ri E[(e^ (x " ) - e ^ (xl) ) 2 ] = -E^^^KX 1 )] 

n 

and we obtain 

OO _2 

hm n £ n = -E[e 2 ^ 7 MP0] E ^^r 1 

= _ e -l e E [e ^(-)] E[e 2MX) 7[v9](x)] 

--exPh/d-e-)*]/^-*]^. 

In order to recognize the limit, let us consider what we have called the white structure 
on N associated with the structure (E, J 7 , (j,, T> , 7) (cf. [5] chapter VI §3, cf. also [7]), 
its square field operator T and its Dirichlet form £ satisfy 

r [ e ij>diV] = _ e 2ij vd N dN 

£[ e *jW] = -lE[e 2l ^ dN fj[ V ] dN] 

which using the Laplace characteristic functional Ee* f f dN = e - ^ 1-6 * ^ dfl , may easily 
be seen to be equal to 



~exp[- J(l-e 2i *)dri j lW*dp 
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what we obtained up to the factor 1/2. In conclusion 



Proposition 15. Let X n be an approximation of X satisfying (H3) on an algebra of 
continuous functions with a local asymptotic Dirichlet form and square field operator^. 
Assuming in addition the weak convergence (X n ,X) =4- (X, X) then the approximation 
N n of the Poisson point process N constructed above satisfies (H3) with the same a n 
and with asymptotic Dirichlet form the so-called white form characterized by its square 
field operator 

r[ f fdN] = J 7 [/]djv. 

II. 9. Stochastic integral. 

We now consider a stochastic integral 

Y= I H s dB s 
Jo 

approximated by the sum 

n— 1 

Y n = V H k (Bk+i - B k ) 

' ^ n n n 

k=0 

(B t ) is a standard Brownian motion defined as the coordinate process of C([0, 1]) 
equipped with the Wiener measure, and H s = H + Jq £ u dB u + r\ u du is an Ito 
process defined on the same space, processes £ and rj are adapted and regular in Malli- 
avin sense. We suppose they satisfy sup i E[|£ t | p + \r]i\ p ] < +oo for some p > 2 and we 
will state their other regularity properties along the calculation. 

In order to obtain the limit expressions we are looking for, we will use several times 
the integration by part formula 

E[uSU] = E[< Du, U > n ] 

(cf. for the notation [5] formula (15) p81). This technique has been already used with 
success by Clement, Kohatsu-Higa and Lamberton [11] to compute, for s.d.e. possibly 
with delay, an estimate of E[<^(Y„) — (p(Y)}, i.e. with our notation, an estimate of 
< A[(p], 1 >. Let us note that this expression which is always equal to < )\[(p], 1 > 
erases the diffusive part of the bias and, since A is here a first order operator, as we will 
see in a moment, this expression writes < A[tp], 1 >=<^[<^], 1] = Ey[Fy/] and reduces, 
when regularity allows an integration by parts, to the form E[G</?]. In the case of Ito 
type s.d.e. under rather general hypotheses, E[y9(F n ) — <^(F)] keeps the same order of 
magnitude for the speed of convergence even when ip is only bounded and measurable 

[2]- 

In this section we attempt to explicit the four bias operators for the above approx- 
imation problem. They occur with the sequence a n = n. 
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a) The local property is satisfied. 

Lemma 4. If £ and r\ satisfy sup t E[|£ t | p + \i]i\ p ] < +00 for some p > 2 then 

nE[\Y n - Y\ 2+a ] ^ Va : 2<2+a<p 

Proof. Let a be s.t. 2 < 2 + a < p, by Burkholder-Davis-Gundy inequality 

E[|y n - F| 2+Q ] < CiE / \H S - H^\ 2+a ds. 

Jo 

Now 

\\H s -H^\\ 2+a < (C7 1 E[/^ 1 eW 5 D' + (E[(/Llk|d«) p ])' 

n n n 

< ( a - M)l (dsii^E^*)? + (a - ^)(su Pt EM^ 
EO^-^n <( s -M)^(C 2 + o(l)). 

Hence nE[\Y n - Y\ 2+a ] < ^(C 3 + o(l)). Q.E.D. o 

It follows that if our test functions algebra V consists of bounded C 2 -functions with 
bounded derivatives, we have lim n nE[|<^(K„) — (p(Y)\ 2+a ] = so that by remark 3 if we 
succeed in proving assumption (H3) with a n = n, the asymptotic Dirichlet form will 
be local. 

b) It follows also if we assume a little bit more for instance that the functions in V are 
C 3 bounded with bounded derivatives, that in the study of A, the expression 

nE[(<p(Y n ) - <p(Y))x(Y)\ 

has the same limit as 

nE[(Y n - Y) V '{Y) X {Y) + \{Y n - Y) V(Y) X (Y)}. 

Similarly, in the study of A, the expression nE[(ip(Y n ) — ip(Y)) 2 } has the same limit as 

nE[(y n -y)y 2 (y)]. 

For simplicity we shall suppose that the functions in V are C°° bounded with 
bounded derivatives. 

c) Study of the symmetric bias operator. 

Let us remark first that in the study of nE[(V r 
and that H be of the form H t = H + f*£ s dB s . 
Indeed, putting K s = J Q S rj u du we have 

nE[( / (K s - K[ns])dB s ) 2 (p' 2 (Y)] < Cn [ E[( / r] u du) 2 ]ds = O(-). 
Jo n Jo J 1 -^ n 



— Y) 2 ip' 2 (Y)} we may suppose r\ = 
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Now for studying nE[(F n — Y) 2 (p' 2 (Y)] we apply Ito's formula to the continuous 
martingale Y n — Y : 

(Y n - Yf = J 1 2 f*(H a - H[ns\)dB s (H t - H M )dB t + J \H S - H^fds 

n n n 

= (1) + (2) 

i) Contribution due to the second term. 
Let us apply once more Ito's formula 

(Jm £,udB u ) 2 = 2 JjL] Jfnt] ^ u dB u ^ t dB t + /[L] 

n n n n 

= (2,1) + (2, 2). 

The contribution of the term (2,2) is nE[J^ J[U £ 2 a dudsip' 2 (Y)] which tends to ^E^ 1 £ 2 dsip' 2 (Y)]. 

n 

The contribution of the term (2,1) is zero. Indeed by integration by parts it is the 
limit of 



2n [ E [ I ^dB^ t D t [ip t2 {Y)\dtds 
Jo Jted jM 



n n 



which by an other integration by part in order to get rid of the stochastic integral, 
gives 

2n [ f f m u D u [i t D t [ V l2 {Y)]]}dudtds 
Jo JhA JM 

n n 

and is 0(^)as soon as E[£ u Z} u [£ t _D t [(/? /2 (F)]]] remains bounded. 

ii) Contribution of the first term. 

We shall show by several integration by parts that this contribution is zero. We are 
concerned by the limit of 

2nE J Q l J \H S - Hjns])dB s (H t - H^)D t [if' 2 (Y)}dt 

n n 

= 2nE /' /'(if, - H^)(H t - H m )D s D t [ip l2 (Y)]dsdt 

n n 

+2nE /o 1 J \H S - H[_ns])D s [H t - H^}D t [ V l2 {Y)}dsdt 

n n 

= (1,1) + (1,2). 
The term (1,1) decomposes in (1,11) + (1,12) with 



(1. 11) = 2n /; /J /^] E[UH t ~ Hw)D u D s D t [ V ' 2 (Y)}}dudsdt 

n n 

(1. 12) = 2n j; jl jU E[i u D u [H t - H^]D s D t y 2 {Y)}]dudsdt 



we have 



l,ll)=2n/ [ [ [ E[i v D v [i u D u D s D t y 2 {Y)}dvdudsdt 
Jo Jo J^ JM 



n n 
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Thus (1,11) is O(^) as soon as the expectation inside is bounded. About (1,12) we get 
still two terms 



(l,12) = 2n/ / / E[U[ D u [Z a ]dB a + Z u l ( M (u))D s D t [ip' 2 (Y)]}dudsdt 

Jo Jo Jhd JM < n V 

n n 

the second one requires — = — and tends to zero. The first one may be written 

2n I I I [ D u [UD a D s D t [ip ,2 {Y)}dadudsdt = 0{-) 
Jo Jo J^ M n 

n n 

Let us come to the term (1,2) which may be written 









(H t - H[nt])[ 


>o Jo 


" Jil 

_ 



dsdt 



still two terms, the second one requires — = — and tends to zero. The first one is 

n n 



2n [ [ E 
Jo Jo 



f S ipdBp f D s [UdB a D t [ip' 2 (Y)} 

[ns] I [nt] 



dsdt 



which may be handled as term (1,1). 

Eventually, supposing £ and n are bounded with bounded Malliavin derivatives up 
to order four, then 

nE[(tp(Y n ) - V (Y)) 2 ] -> \E[j\ldsy'\Y)]. 

d) Study of the theoretical bias A. 

By the remark done in b) above, it is sufficient to study the limit of 

nE[(Y n - Y)ip'(Y) X (Y) + \(Y n - Yf V "(Y) X (Y)]. 

The work is already done for the second term, it suffices to replace (p' 2 (Y) by (p"(Y)x(Y) 
in point c). 

The first term may be written 

nE[fi(H^ - H s )dB sV '(Y)x(Y)} = nEj^H^ - H S )D S [ V > (Y) X (Y)]ds 
= nE f (JiU iudB u + §{ns\ r] u du)D s [(p'(Y)x(Y)]ds 

n n 

= (a) + (b) 

(a) = n JU E[^ u D u D s [ip'(Y)x(Y)\\duds which tends to \ £ E[i s D s D s [ V '(Y) X (Y)]]ds. 
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(&) = f [ S E[r )s D s [ l f'(Y)x(Y)}duds^l f E[rj 9 D 8 [<p'(Y)x(Y)]]ds 
Jo J 1 ^ 1 Jo 

n 

therefore we have 

\im n E[(Y n -Y) V '(Y) X (Y)] 
[ > =\$l E[Z s D s D s [ V '(Y) X (Y)]]ds + \ £ E[ Vs D s [ V '(Y) X (Y)]]ds. 

Hence, with the same hypotheses as for c), the conclusion is : 

\im n nE[(tp(Y n )-tp(Y))x(Y)} 

= \j ® HsDsD^'iYMY)] + r ]s D s y{Y)x{Y)]] ds + \j^ E[^" (Y) X (Y)]ds. 

e) Interpretation of the results. 

If we put ^E[J^ £gds\Y = y] = p(y) the form £[</?] = § / py?' 2 d¥ Y is closable iff the 
measure |pPy has a density satisfying the Hamza condition (cf. [13] pl05). In that 

case the operator A exists and is uniquely defined by 



< A[p], X >= -\E[j\ldsip'{Y)x\Y)\. 



The operator A would be defined by 

< A[<p], X > = lim n nE[(^(r n ) - <p(Y)) X (Y)] 

= \E[j l Q i s D s D s [ V '{Y) X {Y)]ds] + \E[Jlr ]s D s [ V '{Y)x{Y)]ds] 
+\nSo&Wx)'{Y)d8\ - lnio^'(Y)x'(Y)ds] 

provided that the righthand side may be put in the form of the lefthand side. 
Sufficient conditions are easy to be listed, 
(j) By the functional calculus the first term may be written 

l^j\s[(D s [Y]) 2 ^'xr(Y) + D s D s [YW X y(Y)]ds 

it will have the desired form as soon as the measures £ s (D s [Y]) 2 ds.F Y and 
Jq ^ s D s D s [Y]ds.¥ Y will be sufficiently regular to allow an integration by parts, 
(jj) The second term is equal to 

l -E[j\ s dB sV '{Y) X {Y)}. 

(jjj) The third term requires the conditions of an integration by part, 
(jv) The fourth term is < A [ip], x >■ 
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We see that the operator A exists under a quite simple condition involving only 
the law of the pair (J* 1 ^ds, Y). The regularity conditions insuring the existence of the 
operators A or A are more intricate. When they hold, putting 

(12) < A[ V ], X>= l if j\tsD s D s + Vs D s ) [<p>(Y)x(Y)}ds + ±E[ jf * gdaWxYOO] 

we have A = A + A and A = —A + A and by the general theory, hypotheses (HI) to 
(H3) being fulfilled and the Dirichlet form being local, A is a first order operator, as 
may be seen also on the obtained form (12). 

Remark 11. Our approach is direct. But the heaviest part of the proof i.e. the 
proof of 

nEMY n ) - V {Y)f] -> Ie j\ 2 s dsip'\Y\ 
may be considerably shortened if we use a result of weak convergence like 

(V^(Y n - Y), Y) A (-J= J 1 CsdW s , Y) 

with an "extra" Brownian motion W independent of Y and £. 

This gives our results thanks to the uniform integrability of n{Y n — Y) 2 which is 
a consequence of the inequality E[|F n — F| 2+a ] < n~^(Cs + o(l)) established in the 
proof of lemma 4. 

Such a weak convergence result has been obtained long time ago by Rootzen [28 
] for the case where the process H has the form H s = f(B s ,s). This kind of weak 
convergence results for stochastic integrals have been now considerably extended, see 
especially [16], [18], [29], [19], [21], [21], [17], [15]. Our approach to s.d.e. in the next 
section is based on such results. 

11.10. Stochastic differential equations and Euler scheme. 

As we have just explained we will base our approach on results on convergence in 
law, in particular on the article of Jacod and Protter [17]. We consider only the case 
of a continuous semi-martingale in which the main ideas already appear. 

Let X = (X ! ) i=lr ^ be a continuous semi-martingale with values in IR d vanishing 
at zero defined on the stochastic basis (fi,^, (JF 4 ),P). For t G [0, 1] we consider the 
g-dimensional s.d.e. 

(13) dY t = f(Y t )dX t Y = y 

where y G W, f is C 1 from W into R qxd with at most linear growth (|/(a;)| < K(l+\x]) 
denoting |.| the norms on R fc ). It is known that (13) has a unique strong solution. We 
study the resolution of (13) by the Euler scheme : 

dY t n = f(Y^)dX t Y n = y 
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where [nt] is the entire part of nt. 

We denote £7™ = Y t n — Y t the error process. U n as process with values in C([0, 1]) 
tends to zero in probability (as soon as / is locally Lipschitz with at most linear growth 
[17])- 

It is supposed that X = M + A where M is a continuous local martingale vanish- 
ing at zero with values in M, d and A is a continuous finite variation adapted process 
vanishing at zero satisfying 

A\ = J^a l s ds with Jo (a*) 2 ds < +oo a.s. 

< M\Mi > t = f*(? s j ds with f*((%) 2 ds < +oo a.s. 

then for every starting point y and for all function / C 1 with at most linear growth, 
the process converges in law on C([0, 1]) to the solution to 

U?dXi-J2f*(Y t )dZ e t 



dUl 



d q 



3 = 1 k=l 



dp 
dx k 



e=i 



m = o, 



Z being given by 



Z? 



where W is a standard g 2 -dimensional Brownian motion defined on an extension of the 
space independent of X and a is a matrix of processes s.t. (crcr*) 1 - 7 = c y which exists 
as soon as q > d case to which the question may be always reduced. 

The proof consists of the conjonction of theorems 3.3, 5.1 and 5.5 of [17] and their 
proofs. 



In order to study the hypotheses (HI) to (H3) we consider the algebra V of the 
linear conbinations of functions ip defined on C([0, 1]) by 

(fi(Y) = e i<^Y tl >+-+i<u r ,Y tr > U£ eRq Ue [ Q) X | 

and the sequence a n — n. 



l,...,r 



a) Symmetric bias operator. 

We study nE[(y?(y n ) - y?(F)) 2 ]. 



Lemma 5. If for fixed t the sequence n\Y t n — Y t \ 2 = \\/nU^\ 2 is uniformly integrable, 

2' 



(14) 



nE[(<p(Y») - V {Y)f] -> E 



,j=l 1=1 

Proof. Let us argue in the case q = r = 1, the general case being similar. 
nE[(^(y») - V {Y)f] = E[n(Y t n - Y t ) 2 (fi V '(Y t + X(Y t n - Y t ))dXf] 



< E[{n(F« - Y t f - (n(Y t n - Y t f) A a}\\<ff\\ 



/||2 

ooJ 



+E[{(n(r« - Y t f) A a}(fi V '(Y t + \{Y? - Y t ))d\f 
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By the uniform integrability the first term may be made smaller than e > uniformly 
in n by suitable choice of a, then the second term goes to zero by the weak convergence 
of n(Y n — Y) and the convergence in probability of Y t n — Y y to zero. o 



Remark 12. In the classical case of an s.d.e. defining a diffusion process from a 
Brownian motion, if the coefficients are regular, for instance C°° with bounded deriva- 
tives, it is known that v^ll^t" ~~ Yt\\p is bounded for any p G [l,+oo[, the uniform 
integrability of n\Y t a — Y t \ 2 follows. 

Considering that X and W are defined on a product space whose samples are 
denoted u and cD, formula (14) shows that if hypothesis (H3) is verified and if n\U^\ 2 
is uniformly integrable, the limit Dirichlet form satisfies Y t G © and its square field 
operator satisfies 

r\Y?]=&[(u>) 2 ]. 

In other words, the limit process U{uj,Cj) appears to be a gradient in the sense of 
Dirichlet forms of the process Y : we may write 

(15) (Y t f(uj,u) = U t (u J ,u) 

and formula (14) follows by the chain rule. 



The remaining question is whether the form defined on T> by (14) is closable in 
L 2 (C([0, 1]), Py). To this question we have yet only an answer in the simplest case 
where q — 1. When 

dY t = a(Y t , t)dB t + b(Y t , t)dt 
with a, b C 1 with at most linear growth, the process U is given by 

f* a(Y s ,s)a'(Y s ,s) 

u US L dWs 

Jo V2N S 

with t ^ 

N t = exp{ I a' y {Y s , s)dB s - \ I a ,2 (Y s , s)ds + [ b' y (Y s , s)ds}. 
Jo ^ Jo Jo 

Let us denote (£ e ou ,I5> e ou ) the Dirichlet form on the Wiener space of type Ornstein- 
Uhlenbeck with deterministic weight 9, and let us denote D 6 ou its gradient operator 
defined with the auxiliary Hilbert space L 2 ([0, l],dt). We have 

Proposition 16. If the coefficient a satisfies E a'y(Y s , s)ds < +oo and ifa'y(Y s , s) > 
9(s) > 0, hypothesis (H3) is fulfilled. The asymptotic Dirichlet form is the image by Y 
of the form (£ W ,H) W ) defined on the Wiener space by 

D» = {FGC : j\[{Dl[F\{t)f a ^^\dt<^} 
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The proof has been exposed at the Fifth Seminar on Stochastic Analysis, Random 
Fields and Application at Ascona in 2005 and will appear in the proceedings. 
The form (£ W ,U) W ) admits the square field operator 

V w [F] = \j\Dl[F][t)f a ^^dt. 

Putting £ t = ^a'y(Y t ,t) the operator A is given by 

A[p](y) = E[A vl [<p(Y)]\Y = y] 

where ^[(^(F)] = — ^^[l-D^F]], and 5 e ou being the Skorokod stochastic integral 
operator associated with (£^ u ,B>g U ). 

From the concrete point of view of error calculus, the relation 

(16) Y* = U t 

is the most important. It allows to propagate errors by the chain rule and using also, 
in order to manage limit objects, the fact that the operator # is closed what is a con- 
sequence of the closedness of the form. 

b) The theoretical bias operator. 

As in the case of the approximation of a stochastic integral (cf. section II. 9) the 
operator A involves an iterated gradient. 

The main part of the calculation has been performed by Malliavin and Thalmaier 
([24] and [25]) and we adopt their hypotheses : Y is solution of the s.d.e. 

dY t = a(Y t )dB t + b(Y t )dt 

where B is a (d — l)-dimensional Brownian motion and where the matrix a and the 
function b are C°° with bounded derivatives. 

The operator A is given by \im n nE[((p(Y n ) — ip(Y))x(Y)]. Since V consists of 
functions of finite dimensional marginals, we restrict for simplicity to marginals of 
order one and to the case where B and Y are scalar, we have (cf. lemma 4) : 

limnE[Ml7V^))xM]=^ 

n n Z 

P. Malliavin and A. Thalmaier have computed the first term which may be pulled back 
on the Wiener space 

(17) linmE[(>7 - Y t )tp' {Y t ) X {Y t )\ = f E[a 11 (Y s )D s D s F + b x (Y s )D s F + Cl (Y s )F]ds 

n Jo 
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where F = ip f (Y t )x(Y t ) and where an, bi, c\ are functions of the coefficients a, b and of 
their four first derivatives. It should be noted the similarity between the above formula 
(17) and formula (11) obtained for a stochastic integral ((17) reduces to (11) when 
6 = 0). 

The second term is consequence of the preceding results on convergence in law. 

(is) \ E[(y t - - y t )V'(Y t ) x (Y t )] - lnufr"(Y t ) x (Yt)}. 

We see, by formulae (17) and (18) that the operator A is the image by Y of a singular 
distribution operator on the Wiener space. 

We have to conclude that, up to now, the study of the approximation of the solution 
of an s.d.e. by the Euler scheme is far from being achieved : the operator A is yielded 
by the quoted recent works but the existence of the operator A (hence of the Dirichlet 
form) is only shown in a very particular case. 

III. Conclusive comments. 



We focuse in this conclusion on remarks concerning the comparison between deter- 
ministic and stochastic approximation. 

The hypothesis of uniqueness of the approximation of order n. 

Let us consider a situation where given Y the approximation Y n is completely de- 
termined, i.e. for Py-a.e. y, the conditional law of Y n given Y — y, has the form S Vn ( y y 
We call this assumption of uniqueness hypothesis (U). 

Example. Such a hypothesis if often implicitely supposed when numerical results are 
given under the form 

Y 3 = 2.3769 ± 10~ 4 
, v Y 5 = 2.376985 ± 1(T 6 

^ ^ Y 7 = 2.37698534 ± 1(T 8 



and it is underlying the concept of number of significant digits. 

Indeed, let us take the decimal representation of real numbers in [0, 1] : 

oo 

?/ = El^T with anG{0,l,...,9} 

n=0 

If the a n 's are drawn independently uniforly on {0,1,..., 9} the random variable 
Y = Yl^=o KT+ 1 ^ S uniformly ditributed on [0, 1[ and as soon as y is not decimal, 
which is a negligeable set, the expansion of y is unique, so that the above hypothesis 
(U) is fulfilled for the approximation Y n = X]fc=o HF^ • 
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Some martingales satisfy hypothesis (U). Let us consider (fl, A, P) with an increas- 
ing sequence of sub-a-fields B n generated by countable partitions V n of Vl. Then for 
Y e L\ Y n = E[Y\B n ] satisfies (U) with 

Aev n y > 

This happens in particular for Haar systems (cf. [26] chap. Ill §3). 

Proposition 17. Suppose hypothesis (U). If for a n — > +oo and an algebra V, 
a n E[(tp(Y n ) - <p(Y)) X (Y)] ^< A[<p], X >f y ^ eV, V X G L 2 (P y ), 
^/ien (HI) to (H3) /ioW, A = —A =)\ are first order operators and A = 0. 

Proof. The sequence a n (E[<^(y^)|y — y] — <p(y)) is weakly bounded in L 2 (Py) hence 
strongly bounded, i.e. 

a 2 n J (EMY n )\Y = y]- ^y)f F Y (dy) < K. 
Now E[(p(Y n )\Y = y] = <p(r]n(y)), hence 

a n E[( V (Y n ) - V {Y)f] = a n [ (<p( V n(y)) ~ ^(y)) 2 Py(^) < — ^ 0. 

J OL n 

The Dirichlet form is zero, hence it is local and )\ is a first order operator. o 

Examples in part II show that in many probabilistic approximations, hypothesis 
(U) does not hold. The law of Y n given Y = y has a non zero variance. Polya's urn is 
a generic example showing that information at each step cannot be resumed by boxes 
of size ±lCT fc but by standard deviation of laws whose support doesn't go in general 
to zero. 

When we are interested in computing a sample of a random quantity, we have to 
display the result with specifications adapted to the stochastic case. 

The interest of such specifications is particularly clear in infinite dimension when 
we have to compute wbywa path of a process. For example in the GPS or GALILEO 
systems when modelling the ionosphere by a spatio-temporal process, computing a 
sample is necessary to obtain the shift in the signals coming from several satellites. 
The accuracy of this sample is important to get the accuracy of the whole positionning 
system. 

For such numerical computations of sample paths, we suggest that, as much as 
possible, the following specifications be displayed : (i) the law Py of Y, (ii) the sequence 
a n , (iii) the theoretical and practical bias operators A and A. 

Then the operator A, the Dirichlet form and the square field operator follow and 
the approximation Y n (uj) may be (if the Dirichlet form is local) the starting point of 
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an error calculus for the studied model. Non locality of the form, when it happens, is 
also a precious warning to be particularly carefull in the sensitivity analysis. 
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